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Abstract 

Over the past decade, fundamentals of time independent density functional theory for excited 
state have been established. However, construction of the corresponding energy functionals for 
excited states remains a challenging problem. We have developed a method for constructing func- 
tionals for excited states by splitting /c-space according to the occupation of orbitals. In this paper 
we first show the accuracy of kinetic energy functional thus obtained. We then perform a response 
function analysis of the kinetic energy functional proposed by us and show why method of splitting 
the A;-space could be the method of choice for construction of energy functionals for excited states. 
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I. INTRODUCTION 



Time-dependent density functional theory (TDDFT) is now a standard tool for calcu 
lating excited-state energies. Simultaneously parallel efforts have always been made to de 
velop time- independent excited state density functional theory (DFT). These include work 



of Zieglei;^et al. Gumiarsson and Lundqvist 



Pathak [5|, Theophilou jsl, Oliveira, Gross and Kohn |7|, l8|, Nagy |9|, Sen [lO|, Singh and 
Deb Formal development 12-15| akin to ground-state density functional theory, how- 
ever, is more recent. A challenging problem in excited state DFT is the construction of 
excited state energy functional. Without such functional, there is no option but to use the 
existing ground state energy functionals for excited states also. In such a case, the only way 
any information about the excited state can be put into energy is through the density alone 



2[, von Barth [3|, Perdew and Levy 



QB, Nagy S 



However, employing ground state functionals for excited states is both qualitative 



17H19| . Thus there is a need for development of 



incorrect as well as numerically inaccurate 
appropriate excited-state energy functionals. 

Starting point for most of the ground state functionals has been the local-density approx- 
imation (LDA). Building on this, GEA 20|, |2l| and GGA 22|, |23j are further constructed. 
Thus for excited states too one would like to first develop an LDA functional and then im- 
prove upon it by including gradient corrections. In the past, we have proposed a systematic 



method of constructing LDA functionals for excited states. This is done 



/c-space according to the orbital occupation of a given excited state 15 



17 



19 



splitting the 



25j. For 



example, in Figured] we show an excited state where some orbitals (core) are occupied, then 
some orbitals are vacant followed by orbitals that are occupied (shell). The corresponding 
/c-space, also shown in Figured! is constructed according to the orbital occupation such that 
it is occupied from to /ci, vacant from ki to k2 and again occupied from k2 to k^ where ki, 
k2, ^3 are determined by 

kl{r) = 37rVc(r) 
kl{r)-kl{r)=3n'p,{T) 

klir) - ^^^(r) = 37rV.(r) (1) 
in terms of pc, Pv and ps corresponding to the electron densities of core, vacant (unoccupied) 
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(a) 



(b) 



FIG. 1: Orbital occupation of electrons and the corresponding k-space occupation for an excited 
state (Core-Shell) configuration similar to that of homogeneous electron gas (HEG). 



and the shell orbitals. Further, 

i=l 

j=ni+l 
713 

Psir)= Yl kf^"Wr (2) 

i=?i2+l 

where first rii orbitals are occupied, ni + 1 to n2 are vacant followed by occupied orbitals 
from n2 + 1 to n^. The total electron density p(r) is given as 

p(r) = pc(r) + Ps(r) 

or p(r) = pi(r) - p2(r) + psir) (3) 

with pi = Pc,P2 = Pc + Pv and ps = Pc + Pd + Ps- Approximate non-interacting kinetic 
energy functionals for such excited-states are constructed using the ground state energy 
functionals. Thus while up to second order in Vp, the ground state non-interacting kinetic 
energy is approximated as 

i 

where 
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= ^(Svr^)^ / pHr)dr (4) 
is the Thomas- Fermi kinetic energy functional with kpir) = (37r^p(r))^/^ and 

is the second-order gradient correction. The same functional cannot be expected to be as 
good for excited states as it is for the ground states. For excited states, the appropriate 



zeroth order approximation T*(°^ is given 19| by the modified Thomas-Fermi kinetic energy 
functional 

T<'\k,, k,, h] = ^,j (klir) - kl{v) + kl{r)) dr (6) 

where the core term X]r=i is approximated by J kl{r) dr and the shell 

part Er=n2+i {^f^^^l - ki'*'") is approximated by ^ / [A;|(r) - A;|(r)] dr. Similarly, 
the construction is extended to the gradient correction up to the second-order by writing 

72/ p,(r) 72/ pj(r) "'+72/ p^W ^' 



Results obtained with these excited-state functionals for excited-states are far superior [H 
to those calculated with the use of the ground-state functionals . The method is quite general 
and is easily extended to excited state with more than one gap. For example, consider an 
excited state with two gaps with the orbital occupation of electrons shown in Figure [21 
The corresponding split A;-space is occupied from to /ci, from k2 to A;3 and from k^ to k^ 
(core-shell-shell structure). For these excited states, the kinetic energy functional up to the 
second order is 

T<'\ku k,, fcs, fc4, = / (^i(r) - klir) + klir) - kHr) + ^^(r)) dr (8) 

^.(2) ^ 1 / [ iVpiWr _ |Vp2(r)r , |Vp3(r)r _ |Vp4(r)r , iVpsWI^ , ... 
72 7 \ pi(r) p,ir) ^ paW p,ir) ^ psW J ^ ^ 

These functionals for spin-compensated systems (unpolarized) are easily generalized to their 

spin-density counterparts by writing 

^*[pt> Pi] = \iTL^oM + t:,,^oM) (10) 



(a) 



(b) 



FIG. 2: Excited state orbital occupation of electrons for two gap and the corresponding k-space 
occupation similar to that of homogeneous electron gas (HEG). 



where and pi denote the density of the up and the down spin electrons, respectively. 
Results obtained with these functionals for two-gap excited states are shown in Tabled where 
we compare the excited kinetic energies calculated using (|8]) and (Q with the exact Kohn- 
Sham energies obtained by solving the Kohn-Sham equation with Gunnarsson-Lundquist 
parametrization ^ of the LSD for exchange and correlation energy. The density used in 
equation and (11]) is the one obtained from the same Kohn-Sham calculations as done 
for exact kinetic energy. Also shown in Table |T] are the corresponding energies calculated 
from the ground state kinetic energy functional T^°^ and its gradient correction T^'^\ As is 
evident from the table, the ground state based functionals underestimate the exact kinetic 
energies by a large amount whereas the proposed functionals reduce the error significantly. 
This points to the correctness of physics invoked to construct the functionals. 

Similarly, the functional can also be applied to the shell-shell structure where the lowest 
energy orbitals are empty, and the next orbitals are occupied followed by some empty orbitals 
and again orbitals are occupied. The corresponding excited state functional is obtained by 
substituting pi = in equations ([8]) and iQ. This kind of excited states is quite interesting 
because the error of ground state kinetic energy functionals in these is very large. This is 
shown in Table Ull On the other hand excited state functionals of equations ([H]) and ([9]) 
again give energies which are far superior as is evident from Table Ull 

Finally we mention that we have also included the fourth order correction in a manner 
similar to above. The results obtained with the second-order do not change much with the 
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inclusion of the fourth order. 



II. RESPONSE FUNCTION ANALYSIS 

With this accurately constructed functionals based on writing the energy functional E in 
terms of h, A;2, A;3, ■ ■ ■ i.e E[ki, k2, h,---] = E[ki{pi{r)),k2{p2{r)).:h{pz{r))., ■ ■ ■], we ask if it 
is necessary to write the energy functional in terms of pi, P2, Ps or can it be written in terms 
of the density p = (pi — P2 + Ps) directly. Simple answer is that for homogeneous electron gas 
this cannot be done because the density for the ground and any excited-state of the system 
is the same. Thus information about the excited state cannot come from density alone as 
also suggested by the bi-functional nature [3] of excited state functionals. In this paper, we 
demonstrate more rigorously through response function analysis that excited-state energy 
functionals cannot be written in terms of p(r) alone. We show that if we try to write the 
GEA for kinetic energy using total p = (pi — p2 + Ps), it leads to kinetic energy density that 
diverges. We conclude that although it may be possible to write excited energy functional 
in terms of the corresponding density for inhomogeneous systems in some other ways, it is 
most easily done if we instead employ the densities pi, P2, Ps as done above. 

In the following, we first obtain the response function of non-interacting homogeneous 
electron gas and expand it to construct the kinetic energy functional correct up to the second- 
order in Vp(r). We then show that the correction grows exponentially in the asymptotic 
regions of a finite system. 

Suppose one could write the kinetic energy for an excited state of a homogeneous electron 
gas as r*(°)[p]. The system is now made slightly inhomogeneous through a perturbation 



Vext{T^)- The corresponding kinetic energy up to t 



le second order is obtained by expanding 



kinetic energy in density changes and is given as 261-128 



T[p(r)] = T*(°)[p] + IJJ - r'l)Ap(r) Ap(r')rfr dr' (11) 
In the Fourier space, the total energy of the system, using p(q) = f Ap(r)e~'"''''rfr, is 

E[p] = T<%] + I K(q)|p(q)pdq+ ^ J V;.,(q)p(q)rfq (12) 

where K{q) is the Fourier transform of K{\r — r'|). Making the energy stationary with 
respect to p(q) gives p(q) = Vexti^l) K {q)~^ . K{q) is obtained by identifying —K{q)~^ with 
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X^^\<i), where x^°^(q) is the response function of non-interacting homogeneous electron gas 
obtained using perturbation theory and is given as 



k.occ 



k 



k+q 



£0 - £0 
k 



k—q. 



C.C 



where e 



(0) 



and V is the volume over which periodic boundary conditions are applied. 



The response function above for excited state corresponding to Figure [1] is given in a simple 



way as 



XG\ki,k2,k3;q) = Xg' (ki, q) - Xg' (^2] q) + Xg' i^s] q. 



.(0) 



(0), 



(13) 



where Xci^i ^t) is the response function for the ground state with Fermi wave vector k and 
is given as 



xUk;^) 



2 V 2 



2k 



2k 



(14) 



Expanding Xc^iki, k2, k^; q) to order leads to the second order correction for kinetic 



energy. For excited states, small q'-limit gives 



X^'\k,,k2,k3;q) 



ki- k2 + k-s 



SI 
12 



1 



1 



1 

h 



leading to 



-X^°Hq) 



TT 



■ 1 1 1 ■ 

ki k2 ks^ 



iki-k2 + ks) { 12{ki-k2 + k3) 
The excited state kinetic energy functional is then obtained by substituting K[q) in (flTil . 



so that 



r[p(r)]^To[p]4/^ 



|Ap(r 



|Vp(r) 



A;2 + A;3) ' ' 2A J {k^ 
The second term in the expression above is identified as 



k2 + k,y 



■ 1 1 1 

ki k2 ks 



dr{15) 

of functional T*(o) of 



Sp{r)5p{r') 

equation (Q. It is clear from the above that it cannot be written in terms of p = pi— p2 + ps 
alone. To check its correctness if we let ki = k2 and = kp, this term goes over to sp{r)5p(r') 
of functional Ti°^ of equation (jlj) for the ground state. The third term above gives the 
gradient correction to the kinetic energy up to the second order in Vp. Upon substituting 
ki = k2 and k^ = kp this also goes correctly over to ^ J ^^^y^dr, the gradient correction for 
the ground state. However, for excited states, this term is not well behaved for r — t- oo. In 
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this limit, similar to ground state 29l-l3l| densities, excited-state density also varies 32| as 



p(r) ~ exp [— 2-\/— 2 Smax where Smax is the highest occupied orbital energy. Hence k^, ~ 
-2emax)^^'^ 'f] ■ For excited states ki < k2 < so the gradient term in asymptotic 



exp I — o 



limit is proportional to ■'^f^ 



In terms of orbital energies £1,62)^3 corresponding to 



fci, ^2 and /cs, shown in Figure [H this is given by exp [| (V— 2£i — A^J—2e■i) r] . Since, |£i| > 
1^2! > l^^sl , the gradient term diverges in the asymptotic limit for |£:i| > I6I53I if we insist on 
expanding the kinetic energy functional in terms of the density. This diverging behavior of 
the kinetic energy density is also an explicit demonstration of the lack of Hohenberg-Kohn 



theorem for excited states 



Thus additional input is needed |3J] if one wishes to employ 
excited state density for calculation. As shown by our results in the beginning of the paper, 
this is done easily by dealing with pi, P2 and separately right from the beginning. Indeed, 
the exchange energy functional constructed in a similar manner also leads to accurate results 
for excited state energies 



25|. 



III. CONCLUDING REMARKS : 



Exchange and kinetic energy LDA functionals constructed by splitting A;-space for homo- 
geneous electron gas have been shown to be accurate in the past. Analysis of the kinetic 
energy functional based on the response function of an excited-state of HEG suggests that 
it is not possible to construct LDA functionals that are dependent on excited state density 
alone. Splitting the /c-space and working with the corresponding core and shell densities 
separately thus provides a method of constructing density functionals for excited states and 
may pave the way to ground-state like density functional calculations for excited-states too. 
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TABLE I: Approximate kinetic energies (in atomic units) obtained for excited states of some atoms 
through the use of ground state functionals T^^^ and T^^^ + T^"^^ (equations @ and ([5])) and those 
obtained from T*^^^ and 7"*{o)_j_7^*(2) (equations dSD and Q )• Comparison is made with the exact 
kinetic energy resulting from the solution of Kohn-Sham equation for excited states. 

ZEROTH ORDER SECOND ORDER 
Atom 





T{o) 




r-p{Exact) 


T*{0)+j^*{2) 2-{0) +2-(2) 


A^a(ls22s°2/3s°3/) 


131.024 


138.916 


150.803 


152.797 


142.480 


A^a(ls22s°2/3s°3p^) 


133.511 


142.007 


153.377 


155.897 


144.881 




168.132 


179.525 


191.196 


196.374 


181.698 


^/(ls22s°2/3s°3/) 


203.234 


217.535 


231.138 


237.663 


219.206 


5f(ls22s02/3s03/) 


242.410 


260.000 


275.573 


283.740 


260.998 



TABLE IL Approximate kinetic energies (in atomic units) obtained for pure excited states of some 
atoms through the use of ground state functionals T^*^) and T^^^ +T(^) (equations (j3|) and ([5])) and 
those obtained from T*^^'^ and 7"*{o)-)_7"*(2) (equations ([8]) and ([9]) ). Comparison is made with the 
exact kinetic energy resulting from the solution of Kohn-Sham equation for excited states. 

ZEROTH ORDER SECOND ORDER 
Atom 









'-p{Exact) 


2^*(0)_j_2^*(2) 


rpiO) ^ 2^(2) 


B(ls°2s°2p33s°3p2) 


2.794 


6.308 


7.308 


6.571 


3.025 


Af(ls°2s22p33s03p2) 


9.593 


17.393 


20.962 


18.201 


10.345 


F(ls02s°2p43s03p5) 


18.676 


34.908 


36.703 


36.316 


19.788 




31.164 


50.416 


55.958 


52.514 


32.991 


Mg{ls^2s'^2p'^2,s^2>p^) 


50.125 


80.608 


89.788 


84.062 


52.894 


Mg{ls^2s'^2p^2,s^2>p'^) 


61.526 


92.917 


97.778 


96.224 


64.538 


5f(ls°2s22/3s°3/) 


93.290 


139.872 


146.527 


144.983 


97.636 
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